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Abstract 

This paper presents a unified approach for localizing some relevant graph topological in- 
dices via majorization techniques. Through this method, old and new bounds are derived 
and numerical examples are provided, showing how former results in the literature could be 
improved. 
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1 Introduction and preliminaries 

Optimization problems involving Schur-convex or Schur-concave functions have received much at- 
tention in the literature providing some useful applications in different fields: for example they 
became a convenient tool to localize eigenvalues of a real spectrum matrix ([T], |21) ) or more 
generally to obtain bounds for an arbitrary order statistic distribution (jSj). More recently some 
issues related to the structural properties of graphs, characterized in terms of invariants, have been 
explored solving suitable optimization problems via majorization techniques ([3], [7]). 
The aim of this paper is to explore this approach, showing that it is a powerful tool for achieving, 
through a unified scheme, many well-known bounds of some graph topological indices that can 
be expressed in terms of Schur-convex or Schur-concave functions. One major advantage of this 
technique is to set up a strategy for improving these bounds as well as to obtain new ones. We 
start recalling some basic definitions on majorization. More details can be found in |18) . 

Definition 1 Given two vectors y, z G D = {x G M" '■ xi > X2 > ■■■ > Xn}, the majorization order 
y < z means: 

f(y,sk) < (z,sk), fc = l,...,(n-l) 
\(y,s") = (z,s") 

where (•, •) is the inner product in M" and = [1, 1, • • • , 1, 0, 0, • • • 0], j — 1,2, ■ ■ ■ ,n. 

Definition 2 Given a sef 5 C D n {x G M" : (x,s") = a}, a vector x*(S') G S (x^{S) e S) is 
said to be the maximal (minimal) vector in S with respect to the majorization order if ^ 
fx* (5) <x) for each x G 5. 
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Since the upper and lower level sets: 

U{z) = {x e 5a : z < x} , i(z) = {x e 5a : X < z} 

are closed, the existence of maximal and minimal elements of 5 are ensured by its compactness. In 
[T], [5] and [TB], the extremal vectors for some closed subsets of 

Ea = {x e M" ■.xi>X2> ... > > 0} n {x e R" : (x, s") = a} 

have been computed. In particular in [2 the maximal and minimal elements of 

5a = Ea n {x e M" : M, > X, > m„i = 1, • • • ,n} 

where Mi > M2 > • • • > Af„, rrii > 7712, • • • > m„, have been derived, studying also the extremal 
vectors of some subsets of Sa of particular interest: 

5^ = Ea n {x e M" : Af > xi > a;2 > • • • > > m} ( see also [H]) (1) 
5^ = Ea n {x e M" : > a, i = 1, • • • , /i, 1 < h < n, a < ^} { see also [1]) (2) 
5j^ = Ea n {x e M" : a;^ < a, i = /i + 1, • • • ,n, 1 < < (n - 1), a < a}; (see also [Ij) (3) 

5^'') = Ea n {x e R" : Ml > xi > ••• > a;,! > mi, (4) 
M2 > Xh+i > ■ ■ ■ > Xn > m2, rui < Mi, i = 1, 2} . 

Definition 3 A symmetric function (j): A ^ M., A C_ R", is said to be Schur-convex on A (fx < y 
implies 4>{'x.) < 4>{y). IJ in addition 0(x) <4>{y) for x <! y hut x is not a permutation ofy, (j) is said 
to be strictly Schur-convex on A. A function (j) is (strictly) Schur-concave on A if —(j) is (strictly) 
Schur-convex on A. 

Given an interval / C R, and a (strictly) convex function 5 : / — > R, the function 0(x) — X^ILi di^i) 
is (strictly) Schur-convex on /" = / x / x • • • x /. The corresponding result holds if g is (strictly) 

n — times 

concave on 

From the order preserving property of Schur-convex functions, the solution of some constrained 
nonlinear optimization problems of particular interest can be obtained in a straightforward way. 
More precisely, the problem we face is the following: 

max (min) 0(x) 
subject to X G 5 

where 5 C R" is a generic subset which admits maximal vector x*(5) and minimal vector x*(5) 
with respect to the majorization order. If the objective function (j) is Schur-convex, the maximum 
and the minimum are attained in x*(5) and x, (5) respectively; the opposite holds if is a Schur- 
concave function. This enable us to solve problem (P) in a more direct way, avoiding the standard 
approach of Lagrange multipliers. In the next section we illustrate how the procedure above can 
be successfully applied to get upper and lower bounds for some relevant topological indices, which 
found applications in many fields varying from chemistry to network analysis. 
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2 Majorization and bounds for graph topological indices. 



Let us firstly recall some basic graph notations and concepts to be used later. For more details refer 
to [5] and [TU]. Let G = {V,E) be a simple, connected, undirected graph where V — is 

the set of vertices and E QV x V the set of edges. We consider graphs with fixed order |y| = n 
and fixed size = m. Denote by tt = (di, ^2, ••, rfn) the degree sequence of G, where di is the 
degree of vertex Vi , arranged in non increasing order di > d2 > ■ ■ ■ > dn- It is well known that 

n 

= 2m and that if G is a tree, i.e. a connected graph without cycles, to = n — 1. Let A{G) 

1=1 

be the adjacency matrix of G and Ai > A2 > ... > A„ be the set of (real) eigenvalues of A{G). The 
matrix L{G) = D{G) — A{G) is called Laplacian matrix of G, where D{G) is the diagonal matrix 
of vertex degrees. Let ^1 > ^2 — •■• ^ f^n be the eigenvalues of L{G). The following properties of 
spectra of A{G) and L{G) are well known : 

n n 

^A, = 0; ^A^ = 2to; 

i=l 1=1 

E" 2to 
Ml = 2to; > 1 + "1 > ; = 0, fJ-n-l > 0- 
n 

1=1 

Notice that most of the topological indices of graphs are formulated by strictly Schur-convex (Schur- 
concave) functions of the degree sequence as well as the eigenvalues of A{G) or L(G). The corre- 
sponding bounds are generally expressed in terms of size and order of G but they can take also into 
account the degrees of one or more vertices of G. With respect to the degree sequence, one of the 
most popular index is the General Randic index: 

(vi,Vj)^E 

where a is a non zero real number (see [4]). Setting a specific value for a, some very well known 
indices can be obtained: for example, a = 1 corresponds to the Zagreb index Af2(G) (see [19] ) 

while a = — - and a = — 1 to the branching indices (see [10] )• In the recent paper [53], the 

sum-connectivity index has been proposed and in [5] it has been extended to the Generalized sum- 
connectivity index defined as follows: 

xJG)= E ^d, + d,T. 

Note that for a = 1 we obtain the first Zagreb index Mi{G) — X](t)i v )eE (^j + ^i) — X^UiGB "^i • 
Other frequently used indices involve Schur-convex or Schur-concave fimctions of the eigenvalues of 
A{G) and L{G). We recall, among the others: 

1. Energy index: E{G) - Y.l=i \M (0) 

2- sa(G) = Er=/ > « ^ 0, 1 m) 
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3. Kirchhoff index Kf{G) = nJ2"=i — = ns_i{G) ([33], [IS]) 

4. Laplacian Estrada index LEE{G) = ^"^^ e'"- ([H], [IS]) 

In the next section we present theoretical resuhs on upper and lower bounds for the indices previ- 
ously mentioned, recovering well known bounds and improving, in some cases, the existing ones of 
the literature. We postpone to SectionS some numerical examples. 



2.1 General Randic index 

The generalized Randic index can be equivalently expressed as: 

{vi,Vj)eE \{vi,Vj)eE i=l 

Let TT = (rfi, ^2, ■•, rfn) be a fixed degree sequence and x e M™ be the vector whose components are 

n 

df +df, with {vi, Vj) G E. Hence, taking into account that "^df"^^ is a constant, Ra{G) is a Schur 

i=i 

m 

convex function of the degree sequence of G and it is minimal (maximal) if and only /(x) ^ "^xf — 

1=1 

||x||2 is minimal(maximal). Extending the result proved in [Tj,, it is possible to show that 

m n 

^x.= Y id?+d-)=Yd7^" (5) 

i=l (vi,Vj)£E i=l 

and thus X^I^i is a constant. Hence, considering a subset 5* of Sq C M™, where a = X)"=i d'^^^, 
which admits x*(S') and x*{S) as extremal vectors with respect to the majorization order, the 
function / attains its minimum and maximum on S at f{x:t(S)) and f{x*{S)), respectively. The 
general Randic index can be consequently bounded as follows: 

\\Ms)\\i-Edr^' \\x*{s)\\i-j:dr'-' 

< RJG) < . (6) 



Clearly different numerical bounds can be derived, if we characterize suitably the set 5*. The 
structure of S, among the typology of sets given in (Il])-(l4]), will depend on the information available 
on the degree sequence of G. 

For the particular case a — 1, which corresponds to the Zagreb index AI2, this methodology was 
applied in 'T, and, more recently, in [2; where the authors get sharper bounds for the index M2, 
valid for a particular class of graphs having exactly h pendant vertices, i.e. vertices with degree 
one, with degree sequence 

TT = (di, • • • ,dn-h, 1, • • • , !)• (7) 

h 

Starting from (6), the following results allow to compute the Randic index corresponding to a 
generic a, for a class of graphs with h pendant vertices and degree sequence given by (O. Let us 
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assume, furthermore, that the degree sequence tt is such that: 

1 + rf? < + for a > 

dZ-h + dZ-h-i < 1 + for a < ■ 

These conditions, always satisfied for a — —I, assure that the h components of x, corresponding to 
pendant nodes, are separated from the others. Moreover, since the components of x are arranged 
in nonincreasing order, for a > 0, they are assigned to the last h positions of the vector x while, 
for a < 0, to the first h positions of the vector x. Being interested in bounds for the branching 
indices associated to a = — ^ and a = —1, in the sequel we only focus on negative values of a. It 
is easy to see that the choice of a degree sequence equal to JT]), combined with lead us to face 
the following set of type (|4]): 



5 = {x e : 1 + <xh<---<xi< + 1, 

" (Q) 

<+d^<x™<...<a;„+i<d^,, + C-/.-i, (x,s")=^df+i} 

1=1 

whose maximal and minimal elements with respect to the majorization order can be computed by 
corollaries 3 and 10 in |2]. Let us notice that both inequalities in (|6]) are attained if and only if the 
set S reduces to a singleton. This is the case when di = dn-Ui i-e. when all non pendant vertices 
have the same degree (see [2] for some significant examples). 



2.2 Generalized sum-connectivity index 

Let TT be a fixed degree sequence and x £ M™ be the vector whose components are {di + dj), 
{vi,Vj) € E. The function /(x) = X^I^Li strictly Schur-convex for a > 1 or a < 0, while it 

is strictly Schur-concave for < a < 1. Thus, taking into account that J27Li^i ~ X]"=i ^ 
constant, considering a subset S of Ea, where a — J2^=i di, for a > 1 or a < we get 

\\x.iS)C<XaiG)<\\x*iS)C^ (10) 

where ||-||^ stands for the norm. For < a < 1, the bounds are exchanged. 



2.3 Energy index 

Let us point out that E{G) — X]r=i ~ ^"=1 \/^ thus this index is a Schur-concave 
function of the variables Xf, i = 1, • • ■ ,n. Furthermore X]r=i '^i ~ Xi = Xf,i = 1, - ■ ■ ,n. It 

is well known that Ai > If a sharper lower bound for Ai is available, i.e. Ai > k{> ^), we 

have xi > k'^ > (^)^ > The minimal element of the set 

5* = Ea n {x e M" : XI > k^} 



where a — 2m, is given, using Corollary 14 in 2 , by: 



Thus 



E{G) <k + V("-l)(2m-fc2). 
In the case of bipartite graphs, since Ai = — A„, we have xi = xi and thus we face the set 

5 = n {x G R" : X, > i = 1, 2} 

whose minimal element becomes 



2m - 2fc2 



2m - 2fc2 



n-1 



and as a consequence 



£;(G) < 2fc + v/(ri - 2) (2m, - 2fc2). 

For fc = ^ we get the bounds in [n] and [12], while for fc = the bounds in |27]. Different 

choices of fc allow to obtain all the other bounds in [16] and a new bound for E{G) as soon as a 
sharper lower bound of Ai is available. 



2.4 Laplacian indices 

1. Let Sq(G) — 'Y^i^x Mf jQ! ^ 0, 1 be the index given by the sum of the a-th power of the non 
zero Laplacian eigenvalues and consider the set 

S" = Sa n {/x e R"^i : Ml > 1 + di} 

where a = Y^lzl fi^ = 2m and Sa C R"-i. Since < (i + rf^)^ by Corollary 14 in ^ it 
follows that the minimal element of S is 



l + di, 



2m — 1 — di 



2m — 1 — di 



Taking into account the Schur-convexity or Schur-concavity of the functions Sa (G) , the bounds 
in |26j . Theorem 3 can be easily derived. The same approach can be used to find the bounds 

for bipartite graphs given in [26], Theorem 5, observing that in this case /Lt^ > 2\J'^^^^. 

The bounds in [5S] , Theorem 3 can be also improved taking into account further information 
on the localization of the eigenvalues. For instance, since ^2 — '^2 (see [5;)j we can consider 
the set 

S = T,an{fj,eR"-^ : fj.x>l + di, fi2>d2} (11) 
whose minimal element, computed applying Theorem 8 in [2], is given by 



x^,{S) 



1 + di,d2, 



2m — 1 — di — ^2 



2m — 1 — di — ^2 



n — 3 
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Thus 



^Sa(G) < + + 0<a<l 



(12) 



1 



2. Notice that bounds on the KirchhofF index Kf{G) — nJ2i=i — — ns-i{G) can be easily 

Mi 

derived by bounds on s_i(G'). In particular, using the set PT|) we get 



KfiG) > n 



1 



1 



(n-3)2 



1 + di d2 {2m - 1 ~ di ~ d2) 



(13) 



3. The Laplacian Estrada index LEE{G) = X]r=i ^ Schur-convex function of the Laplacian 
eigenvalues. Thus, proceeding as before we can easily recover the bounds in ^28^, Remark 2. 
Furthermore, if we use the set ([TH) . we find the bound: 



LEE{G) > 1 + e^+'^i + e''^ + {n - 3)e 

3 Numerical examples 



2m-l-di -do 



(14) 



In this section we provide some numerical examples to illustrate our results and to show that, at 
least for the considered cases, we achieve sharper bounds with respect to recent literature. 

i) General Randic index 

Let us consider a tree T with the degree sequence tt = (5, 3, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1). We 
start exploring the Randic index i?(T'), with a — — — . The comparison (see Table 1) with bounds 
in fT3| shows that bounds obtained by applying Corollary 3 and 10 in [2] always perform better: 



Ref. 


Lower 


Upper 


m 


5.891 249 (see Theorem 2.6) 


7.486 884 (see Theorem 2.24) 


m 


6.768 317 


7.064497 



Table 1: Lower and upper bounds for R{T). 
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For the case with a = — 1 and the tree T with the same degree sequence, making a comparison 
with the result in [13_, Theorem 3.7, we get 



Ref 


Lower 


Upper 


13 


1 


4.888889 


m 


3.2 


3.666667 



Table 2: Lower and upper bounds for i?_i(r). 

ii) Generalized sum-connectivity index 

Let us consider a graph with degree sequence n = (3, 2, 2, 1) and pick a = —5. Since, in this case, 
n is even, we make a comparison (see Table 3) between the results in [B] and those obtained from 
Corollary 3 in [2]. We can present here also the case when n is odd, always choosing a = —5. For 



Ref. 


Upper 


m 


9.207015 X 10"^ (see Theorem 1) 


m 


5.075226 X 10"^ 



Table 3: Upper bounds for Xai^) with n even. 



a tree with the following degree sequence n — (4, 4, 4, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1), computing the bounds 
(see Table 4) we have: 



Ref. 


Upper 


i 


2.487446 x 10"^ (see Theorem 1) 


m 


4.94124 x 10"^ (see Corollary 3) 



Table 4: Upper bounds for Xai^) with n odd. 



iii) Laplacian indices 

Let us consider a graph whose degree sequence is tt = (3,2,2,2,1). 
a) Kirchhoff index 

For this index, making a comparison (see Table 5) between the result obtained by applying the 
lower bound obtained in equation (fTS]) and the one derived from [25], Proposition 1, we have 



Ref 


Lower 


Eg 


12.5 


formula (1131) 


13.75 



Table 5: Lower bounds for Kf{G). 
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b) Laplacian Estrada index 



In this case, we draw a comparison between the application of formula (|14p and the result 
provided in [32], Proposition 3.4 (see Table 6) 



Ref 


Lower 


[22] 


59.90761 


formula (fTH) 


68.42377 



Table 6: Lower bounds for LEE{G). 



4 Conclusion 

In this paper we present a unified approach for localizing some relevant graph topological indices, 
based on the optimization of Schur-convex or Schur-concave functions. Our results have been 
derived through the characterization of extremal vectors with respect to the majorization order, 
under suitable constraints. We have shown that classical results can be recovered and sometimes 
improved. Our theoretical approach paves the way to find new bounds, by taking advantage of 
further information that can be extracted from the graph G, that allows us to define tighter sets 

5 for performing the optimization problem. Finally, other topological indices can be localized 
whenever they can be expressed as Schur-convex or Schur-concave functions. 
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